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We consider a parametrically driven damped discrete nonlinear Schrodinger 
(PDDNLS) equation. Analytical and numerical calculations are performed to 
determine the existence and stability of fundamental discrete bright solitons. 
We show that there are two types of onsite discrete soliton, namely onsite 
type I and II. We also show that there are four types of intersite discrete 
soliton, called intersite type I, II, III, and IV, where the last two types are 
essentially the same, due to symmetry. Onsite and intersite type I solitons, 
which can be unstable in the case of no dissipation, are found to be stabilized 
by the damping, whereas the other types are always unstable. Our further 
analysis demonstrates that saddle-node and pitchfork (symmetry-breaking) 
bifurcations can occur. More interestingly, the onsite type I, intersite type I, 
and intersite type III-IV admit Hopf bifurcations from which emerge peri- 
odic solitons (limit cycles) . The continuation of the limit cycles as well as the 
stability of the periodic solitons are computed through the numerical contin- 
uation software Matcont. We observe subcritical Hopf bifurcations along the 
existence curve of the onsite type I and intersite type III-IV. Along the exis- 
tence curve of the intersite type I we observe both supercritical and subcritical 
Hopf bifurcations. 



1 Introduction 

In this paper, we consider a lattice model governed by a parametrically driven 
damped discrete nonlinear Schrodinger (PDDNLS) equation 

i(j>n = -eA 2 <p n + A<p n + 7</>„ - io>4> n - a\4> n \ 2 4> n . (1) 

In the above equation, cj> n = <p n (t) is a complex- valued wave function at 
site n, the overdot and the overline indicate, respectively, the time derivative 
and complex conjugation, e represents the coupling constant between two 
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adjacent sites, Z\ 2 0„ — ^n+i — 2<^> n +<^ n _i is the one-dimensional (ID) discrete 
Laplacian, 7 is the parametric driving coefficient with frequency A, a is the 
damping constant, and a is the nonlinearity coefficient. Here we confine our 
study to the case of focusing nonlinearity, i.e., by setting positive valued a 
which then can be scaled, without loss of generality, to a = +1. 

In the absence of parametric driving and damping, i.e., for 7 = and 
a = 0, Eq. (1) reduces to the standard discrete nonlinear Schrodinger (DNLS) 
equation which appears in a wide range of important applications [1]. It is 
known that the DNLS equation admits bright and dark solitons with focusing 
and dcfocusing nonlinearities, respectively. The stability of discrete bright 
solitons in the DNLS system has been discussed, e.g., in Refs. [3, 4, 5], where 
it was shown that one-excited-site (onsite) solitons are stable and two-excited- 
site (intersite) solitons are unstable, for any coupling constant e. Moreover, the 
discrete dark solitons in such a system have also been examined [6, 7, 8, 9, 10]; 
it is known that intersite dark solitons are always unstable, for any e, and 
onsite solitons are stable only in a small window of coupling constant e. 

Furthermore, the parametrically driven discrete nonlinear Schrodinger 
(PDNLS) equation, i.e, Eq. (1) with no damping (a = 0), has been stud- 
ied in [11] for the case of focusing nonlinearity, where it was reported that an 
onsite bright discrete soliton can be destabilized by a parametric driving. The 
study of the same equation was extended for the other variants of discrete 
solitons in [12], showing that a parametric driving can not only destabilize 
onsite bright solitons, but also stabilize intersite bright discrete solitons as 
well as onsite and intersite dark discrete solitons. In the latter, the PDNLS 
model was particularly derived, using a multiscale expansion reduction, from a 
parametrically driven Klein- Gordon system describing coupled arrays of non- 
linear resonators in micro- and nano-electromechanical systems (MEMS and 
NEMS). 

The discrete nonlinear Schrodinger equation with the inclusion of para- 
metric driving and damping terms as written in Eq. (1) was studied for the 
first time, to the best of our knowledge, by Hcnnig [22] focusing on the exis- 
tence and stability of localized solutions using a nonlinear map approach. He 
demonstrated that, depending upon the strength of the parametric driving, 
various types of localized lattice states emerge from the model, namely peri- 
odic, quasiperiodic, and chaotic breathers. The impact of damping constant 
and driving (but external) in the integrable version of the DNLS system, i.e., 
the discrete Ablowitz-Ladik equation, has also been studied [23] which con- 
firmed the existence of breathers and multibrcathcrs. In deriving Eq. (1), one 
can follow, e.g., the method of reduction performed in [12] by including a 
damping term in the MEMS and NEMS resonators model. 

On the other hand, the continuous version of the PDDNLS (1), i.e., when 
4> n ~ 4> an d — £^2^n ~ dx4>i was numerically discussed earlier in [24] result- 
ing in a single-soliton attractor chart on the (7, a)-plane from which one may 
determine the regions of existence of stable stationary solitons as well as sta- 
ble time-periodic solitons (with period- 1 and higher). Instead of using direct 
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numerical integration as performed in the latter reference, Barashcnkov and 
co-workers recently proposed obtaining the time-periodic one-soliton [25] and 
two-soliton [26] solutions as solutions of a two-dimensional boundary-value 
problem. 

Our objective in the present paper is to examine the existence and sta- 
bility of the fundamental onsite and intersite excitations of bright solitons in 
the focusing PDDNLS (1). The analysis of this model is performed through 
a perturbation theory for small e which is then corroborated by numerical 
calculations. Such analysis is based on the concept of the so-called anticon- 
tinuum (AC) limit approach which was introduced initially by MacKay and 
Aubry [27]. In this approach, the trivial localized solutions in the uncoupled 
limit e = are continued for weak coupling constant. Moreover, our study here 
is also devoted to exploring the relevant bifurcations which occur in both sta- 
tionary onsite and intersite discrete solitons, including time-periodic solitons 
emerging from Hopf bifurcations. For the latter scheme, we employ the nu- 
merical continuation software Matcont to path-follow limit cycles bifurcating 
from the Hopf points. 

The presentation of this paper is organized as follows. In Sec. 2, we firstly 
present our analytical setup for the considered model. In Sec. 3, we perform 
the existence and stability analysis of the discrete solitons through a pertur- 
bation method. Next, in Sec. 4, we compare our analytical results with the 
corresponding numerical calculations and discuss bifurcations experienced by 
the fundamental solitons. The time-periodic solitons appearing from the Hopf 
bifurcation points of the corresponding stationary solitons are furthermore 
investigated in Sec. 5. Finally, we conclude our results in Sec. 6. 

2 Analytical formulation 

Static localized solutions of the focusing system (1) in the form of <f>„ — 
u n , where u n is complex valued and time- independent, satisfy the stationary 
equation 

-eA 2 u n + Au n + 7«„ - iau n - \u n \ 2 u n = 0, (2) 

with spatial localization condition u n — > as n —> ±oo. We should notice that 
Eq. (2) (and accordingly Eq. (1)) admits the reflection symmetry under the 
transformation 

u n -> -u n . (3) 

Following [24, 25, 26], we assume that both the damping coefficient a and 
the driving strength 7 are positive. For the coupling constant e, we also set 
it to be positive (the case e < can be obtained accordingly by the so-called 
staggering transformation u n — > (— l) n u n and A — > (A — 4e)). The range of 
the parameter A is left to be determined later in the following discussion. 

In the undriven and undamped cases, the localized solutions of Eq. (2) can 
be chosen, without lack of generality, to be real- valued (with A > 0) [3] . This 
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is no longer the case for non-zero 7 and a in the stationary PDDNLS (2), 
therefore we should always take into account complex- valued u n . By writing 
u n = a n + ib n , where a n ,b n € R, and decomposing the equation into real and 
imaginary parts, we obtain from Eq. (2) the following system of equations: 

-sA 2 a n + (A + j)a n + ab n - (a 2 + b 2 n )a n = 0, (4a) 

-eA 2 b n + {A- j)b n - aa n - (a 2 n + b 2 n )b n = 0. (4b) 

Thus, the solutions of Eq. (2) can be sought through solving the above system 
for a n and b n . 

Next, to examine the stability of the obtained solutions, let us introduce 
the linearisation ansatz </> n = u n + Se n , where S < 1. Substituting this ansatz 
into Eq. (2) yields the following linearized equation at 0{5): 

ie n = -sA 2 e n + Ae n + je n - iae n - 2\u n \ 2 e n - u 2 e n . (5) 

By writing e„ = ?7„e 1 "* + C„e _4 "*, Eq. (5) can be transformed into the eigen- 
value problem (EVP) 

' eA 2 -A + ia + 2\u n \ 2 ul~1 

l-u 2 n -eA 2 +A-ia-2\u n \ 2 

The stability of the solution u n is then determined by the eigenvalues ui, i.e., 
u n is stable only when Im(w) > for all eigenvalues co. 

As the EVP (6) is linear, we can eliminate one of the eigenvectors, for 
instance £„, so that we obtain the simplified form 

[r+(e)r_(e) - 4(a„6„) 2 ] r, n = {u- ia) 2 Vn , (7) 

where the operators C + (e) and C-(s) are given by 

£+(e) = -eA 2 - (a 2 n + 3b 2 n - A + 7), 
J C_(e) = -eA 2 -(3a 2 n + b 2 n -A- 1 ). 





Vn 


= LU 


Vn 









■ (6) 



3 Perturbation analysis 

Solutions of Eq. (2) for small coupling constant e can be calculated analytically 
through a perturbative analysis, i.e., by expanding u n in powers of e as 

u n ^u^+eu^+s 2 ui? + -- - . (8) 

Solutions u n — Un ^ correspond to the case of the uncoupled limit e = 0. 
For this case, Eq. (2) permits the exact solutions Un^ = al°' + ibn^ in which 
can take one of the following values 
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(0,0), (sA + ,-sB_), (sA_,-sB + ), (9) 



where 






~a 2 )(A± ^ 7 2 


-a 2 ) 




'(7± V7 2 


-a 2 )(/l=F yV 


-a 2 ) 



27 

and s = ±1. Due to the reflection symmetry (3), we are allowed to restrict 
consideration to the case s = +1. 

Following the assumption j,a > 0, we can easily confirm that nonzero 
(A+,— B-) and (A-,— B + ) are together defined in the following range of 
parameters 

yl>7>a>0. (10) 

In particular, when -f = a, the values of (A + , —B_) are exactly the same as 
(A_,-B + ). 

Once a configuration for w„ is determined, its continuation for small e can 
be sought by substituting expansion (8) into Eq. (2). In this paper, we only 
focus on two fundamental localized solutions, i.e., one-excited site (onsite) and 
in-phase two-excited site (intersite) bright solitons. Out-of-phase two-excited 
site modes also referred to as twisted discrete solitons (see, e.g., [28]), which 
exist in the model considered herein, are left as a topic of future research. 

Next, to study the stability of the solitons, we also expand the eigenvector 
having component r\ n and the eigenvalue to in powers of e as 

Vn =riW +er]W +0(e 2 ), lj = o/°> + euo^ + 0(e 2 ). (11) 

Substituting these expansions into Eq. (7) and collecting coefficients at suc- 
cessive powers of e yield the 0(1) and 0(e) equations which are respectively 
given by 

= 0, (12) 
LriP = h, (13) 

where 

L = £+(0)£_(0) - 4( a (°)&(°)) 2 - (a/°> - t a)\ (14) 
/„ = [C-(0)(A 2 + 2a™ a™ + 6&<?>6£>) + C + (0)(A 2 + 2b^b^ + 6a^a^) 

+8a(°)&(%(°>&« + a^) + 2cA« 2ia^] rtf* . (15) 

One can check that the operator L is self-adjoint and thus the eigenvector 
h = col(..., Vn-i' 7 !™^ '^n+i' •••) i s m tne null-space of the adjoint of L. 
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From Eq. (12), we obtain that the eigenvalues in the uncoupled limit e = 

are 

c4 0) = ±V A 2 - 7 2 + ia, (16) 

and 

4 0) = ±/c+(0)£_(0) - 4(al 0) 6l 0) ) 2 + ia, (17) 

which correspond, respectively, to the solutions ul ' = (for all n) and 
= al ' + ib^ 7^ (for all n) . For bright soliton solutions having boundary 
condition u n — > as n — > ±oo, the eigenvalues and cj^ have, respec- 
tively, finite and infinite multiplicities which then generate a corresponding 
discrete and continuous spectrum as e is turned on. 

Let us first investigate the significance of the continuous spectrum. By 
introducing a plane- wave expansion r\ n = fie lKn + ve~ %Kn , one can obtain the 
dispersion relation 

u = ±v/(2£(cosk - 1) - A) 2 - 7 2 + ia, (18) 
from which we conclude that the continuous band lies between 

ujl — ±V A 2 - 7 2 + ia, when k = 0, (19) 

and 

Lou = ±\/A 2 -7 2 + 8e(A + 2e) + ia, when k = n, (20) 

From the condition (10), one can check that all the eigenvalues u> £ ±[u>l,Wu] 
always lie on the axis Im(cj) = a > for all e, which means that the continuous 
spectrum does not give contribution to the instability of the soliton. Therefore, 
the analysis of stability is only devoted to the discrete eigenvalues. Discrete 
eigenvalues that potentially lead to instability are also referred to as critical 
eigenvalues. 

3.1 Onsite bright solitons 

When e = 0, the configuration of an onsite bright soliton is of the form 

=0 for n 0, =A + iB, (21) 

where (A, B) ^ (0,0). From the combination of nonzero solutions (9), we can 
classify the onsite bright solitons, indicated by the different values of (A, B), 
as follows: 



(i) Type I, which has (A, B) = (A+, -BJ), 

(ii) Type II, which has (A, B) = (A_, -B+), 
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and 



which we denote hereinafter by u„/±} and u n ^ T y, respectively. 

The continuation of the above solutions for small e can be calculated from 
the expansion (8), from which one can show that an onsite soliton type I and 
type II, up to 0(e 2 ) } are respectively given by 

; ' ^^7& n = ~^ (22) 

0, otherwise, 

'''<=> = t£^& n =~^ (23) 

0, otherwise. 

In particular, when a = 7, the onsite type I and type II become exactly the 
same. 

To examine the stability of the solitons, we need to calculate the corre- 
sponding discrete eigenvalues for each of type I and type II, which we elaborate 
successively. 

3.1.1 Onsite type I 

One can show from Eq. (12) that at e = 0, an onsite bright soliton type I has 
a leading-order discrete eigenvalue which comes as the pair 

= ±VV + ia, (24) 

where 

V = AA^/j 2 -a 2 + 4 7 2 - 5a 2 . (25) 

The eigenvector corresponding to the above eigenvalue has components rjn ^ = 
for n 7^ and r/^ = 1. 

We notice that V can be either positive or negative depending on whether 
a ;=s ath, where 



ath 



= %\lh~i 2 - 2A 2 + A\J A A 2 + f^. (26) 



Therefore, the eigenvalue wj^j can be either 



4±} 

for the case a < a t h, or 



= ±\f AA^f^a 2 + Aj 2 - 5a 2 + ia, (27) 



u 2 fl } =i(a± \lha 2 -4/1V7 2 -a 2 - 4 7 2 ^ , (28) 
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for the case a t h < a < j. 

The continuation of the eigenvalues (29) and (30) for nonzero e can be 
evaluated from Eq. (13) by applying a Fredholm solvability condition. As the 
corresponding eigenvector has zero components except at site n = 0, we only 
need to require /o = 0, from which we obtain the discrete eigenvalue of w n {±} 
for small e, up to 0(e 2 ), as follows. 

(i) For the case a < a t h- 



+ia. 



u; 1{±} = ±\/AA^ 1 2 -a 2 + A 1 2 -5a 2 ± (4 ^ 7 ' a2)S 

J A A ^7 2 - a 2 + 47 2 - 5a 2 

(29) 

(ii) For the case a t h < a < 7: 



w 2{±} 



= i [a ± y ha 1 - AA^J ^ 2 - a 2 - A<y 2 =p 



(4 \Jl 2 - ol 2 )e 



\j 5a 2 - 4/1 yV - a 2 - 4 7 2 
(30) 

We should note here that the above expansions remain valid if ±V are 0(1). 

Let us now investigate the behavior of the above eigenvalue in each case. In 
case (i), the imaginary part of ^i{±} (i-e., when e = 0) is a, which is positive. 

We also note that l^i^l ^ when a ^ a cp , where 

1 



\/257 2 -/l 2 . (31) 

As e increases, the value of |wi{±} | also increases. As a result, the eigenvalues 
w i{±} wm collide either with the upper band {oju) of the continuous spec- 
trum for a < a cp , or with the lower band (wl) for a cp < a < ath- These 
collisions then create a corresponding pair of eigenvalues bifurcating from the 
axis Im(wi{ ± i.) = a. This collision, however, does not immediately lead to 
the instability of the soliton as it does for a = [11, 12]. In addition, the 
distance between wi^j and oj^ increases as a tends to 0, which means that 
the corresponding collisions for smaller a happen at larger e. From the above 
analysis we hence argue that for a < ath and for relatively small e, the onsite 
soliton type I is stable. 

In case (ii), it is clear that \J ha 2 — AA^J^ 2 — a 2 — Aj 2 < a which implies 

< min(Im(w2{2})) ^ a > ^ e fitter indicates the soliton is stable at e = 
0. As e increases, both max(Iin(w 2 {±))) and min(Im(a;2{±})) tend to a at 
which they finally collide. From this fact, we conclude that for small e and for 
c*th < ct < 7, the soliton remains stable. In particular, when a — 7, we have 
min(Im(w 2 {±})) = for all e, which then implies that the soliton is always 
stable. 
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3.1.2 Onsite type II 

Performing the calculations as before, we obtain that the discrete eigenvalue 
(in pairs) of an onsite bright soliton type II is given, up to 0{e 2 ) 1 by 



a ± \J 4^1/7^"- a 2 - 47 2 + 5a 2 ± 



(4V7 2 a2 ) £ 



\J 4:A s /^f r - a 2 - 47 2 + 5a 2 / 



(32) 

Again, we should assume that the term (4/1 y^ 2 — a 2 — A"/ 2 + 5a 2 ) in the 
above expansion is 0(1). 

When a < 7, we simply have \J AAy/j 2 — a 2 — 4"/ 2 + ha 2 > a, from which 

we deduce min(Im(u;| < ^|)) < 0, meaning that at e — the soliton is unstable. 
In fact, as e increases, the value of m.m(lm(uj^ T y)) decreases. Therefore, in 
this case we infer that the soliton is unstable for all e. 

When a = 7, by contrast, the value of min(Im(a;{ =F j)) is zero for all e, 
which indicates that the soliton is always stable. In fact, the stability of an 
onsite type II in this case is exactly the same as in type I. This is understand- 
able as the onsite type I and type II possess the same profile when a = 7. 



3.2 Intersite bright solitons 

Another natural fundamental solution to be studied is a two-excited site (in- 
tersite) bright soliton whose mode structure in the uncoupled limit is of the 
form 

( A + iB a , n = 0, 
= I A 1 +iB 1 ,n = l, (33) 
[ 0, otherwise, 

where (Ao,Bo) 7^ (0,0) and (A\,B\) 7^ (0.0). The combination of the nonzero 
solutions (9) gives the classification for the intersite bright solitons, indicated 
by different values of (A , B ) and (Ai,Bi), as follows: 

(i) Type I, which has {A ,B ) - (Ai.Bi) = {A+, 

(ii) Type II, which has (A , B ) = (A^Bx) = (A_,-B + ), 

(iii) Type III, which has (A Q ,B Q ) = (A + ,-B_) and (A 1 ,B 1 ) =(A_,-B + ), 

(iv) Type IV, which has (A ,B ) = (A_,-B + ) and {A 1 ,B 1 ) = (A + , -B_). 

Let us henceforth denote the respective types by u n ^±±y, u n { TT y, u n {± T y, 
and u n{T±} . 

From the expansion (8), we obtain the continuation of each type of solution 
for small e, which are given, up to order e 2 , by 



10 



M. Syafwan, H. Susanto, and S. M. Cox 



M±±} 



{{A 

(A 



iB- 
iB. 

( L A + -iB-)e 



1 (A+- 

1 (A+-iB-)s 

2 A+y/- 



iB-)e n 

' , n = 0, 
n = 1, 



/7* 



= < 



iB 4 



0, 
(^- 

(A_-i_B + )e 
A—\f' 

0, 



1 (A_-iB + )e 

2 yl_^/ 7 2_„2 ' 



n = -l,2, 
otherwise, 

n = 0, 



LO +) + 2 . / .,2 -2 ' " ~~ ij 



A — \f 



7 , 



7" 



M=F±} 



(A_ - iB + ) + i 

(A_-iB_)e 
A—<\/-y 2 —a 2 ' 

o, 

(A_-zB_)e 
A—^J 7 2 — a 2 ' 

(A_ - zB+) + | 



(■A+-iB-)e 



(■4--iB+) e 



7(yl- A /7 2 -Q 2 ) 



(■4_-»g- t 



2 7( J 4-V7 2_ ^ r ) 
^ x _|_ 1 (A+-iB-)e 

B+)e 



7(A+V7 2 -" 2 ) 



o, 



n = -l,2, 
otherwise, 

n = —1, 
, n = 0, 
, n = 1, 

n = 2, 

otherwise, 

n = —1, 
, n = 0, 
, n = 1, 

n = 2, 

otherwise, 



where 



-4± = 2jA± + {A± v / - f 2 -a 2 )A T , 
B± = 2 7 S T -(A± vV - a 2 )B±. 



(34) 



(35) 



(36) 



(37) 



(38) 
(39) 



All solutions above are defined on the region (10) and exhibit the same profiles 
when a = 7. One can check that intersite type III and IV are symmetric, thus 
they should really be considered as one solution. However, we write them here 
as two 'different' solutions because, as shown later in the next section, they 
form two different branches in a pitchfork bifurcation (together with intersite 
type I). 

Let us now examine the stability of each solution by investigating their 
corresponding discrete eigenvalues. 



3.2.1 Intersite type I 



By considering Eq. (12) and carrying out the same analysis as in onsite type 
I, we obtain that the intersite type I has the double leading-order discrete 
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eigenvalue 



(o) 



Wl {±±} = ± V 4yl V7 2 - a 2 + 47 2 - 5a 2 + ia, (40) 
for a < a t h, and 

w 2 f2 ±} = i (a ± ^a 2 -4/1^72- a 2 - 4 7 2 ^ , (41) 

for a t h < ol < 7- The corresponding eigenvector for the above eigenvalues has 
components 77!°^ = for n / 0, 1, ?7q ^ 7^ 0, and 77^°' / 0. 

One can check, as in onsite type I, that the position of wij^j., relative to 

Wp' depends on whether a a cp = g \/257 2 — /l 2 , i.e., the value of 

is greater (less) than | when a is less (greater) than a cp . 

The next correction for the discrete eigenvalues of an intersite type II can 
be calculated from Eq. (13), for which we need a solvability condition. Due to 
the presence of two non-zero components of the corresponding eigenvector at 
n = 0,1, we only require /o = and /1 =0. Our simple analysis then shows 
t/q - 1 = ir^ -* from which we obtain that each of double eigenvalues (40) and 
(41) bifurcates into two distinct eigenvalues, which are given, up to order e 2 , 
as follows. 

(i) For the case a < a t h- 

wii {±±} = ±V / 47l v /r^ 2 " + 47 2 ~5a 2 ± (2V7 2 - a 2 )e + .^ 

J 4A^-f 2 - a 2 + 47 2 - 5a 2 

(42) 

lo 12{±±} = ±y / 4/ly/ 7 2 - a 2 + 4 7 2 - 5« 2 T 1 2 ( ^+V7 2 -^ 2 )£ + .^ 

y 4/1 Vt 2- a 2 + 47 2 - 5a 2 

(43) 

(ii) For the case a t h < a < 7: 



w 2i{±±} = i I a ± y 5a 2 — 4/1 y/j 2 — a 2 — 4"/ 2 =p 



± ^ 5a 2 - 4Ay/j 2 - a 2 - Aj 2 ± 



(2v^ 



w 22{±±} - * 



y/ 5a 2 - 4/1 A/7 2 - a 2 ' - 4~f 2 
(44) 

2(/l+ v / 7 2 - Q; 2 )g 

y^ 5a 2 - 4/1 \/7 2 - a 2 - 47 2 
(45) 



As before, we assume here that the terms ±(AAy / -f 2 — a 2 — 4-f 2 + 5a 2 ) are 
O(l) so that the above expansions remain valid. 

Let us first observe the behavior of the eigenvalues in case (i). In the 
uncoupled limit e = 0, the imaginary part of wn|^ ± j = U)i2^±±y is a > 
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which indicates that the soliton is initially stable. When e is turned on, 
the value of increases but |wi2{±±}| decreases. Therefore, we can 

determine the mechanism of collision of these two eigenvalues with the inner 
or outer boundary of continuous spectrum (wj, or wu) as follows. 

• For a < a cp , the first collision is between wi 2 {±±} and cou- Because ojjj 
moves faster (as e is varied) than uu{±±y, the next collision is between 
these two aforementioned eigenvalues. 

• For a > a cp , the mechanism of collision can be either between uji 2 {±±y 
and Wi, or between wi2{±±} and itself. 

All of the mechanisms of collision above generate new corresponding pairs 
of eigenvalues bifurcating from their original imaginary parts, which is a. 
Yet these collisions do not immediately cause an instability, because a > 0. 
Therefore, we may conclude that for sufficiently small e and for a < ath, an 
intersite bright soliton type I is stable. 

Next, we describe the analysis for the eigenvalues in case (ii). When 

e = 0, we have < min(Im(w2i /2 ± i )) = mm (I m ( w 22{2 ± })) < a. As e is 
increased, min(Im(w2i{±±})) increases but min(Im(u;22{±±})) decreases. The 
latter then becomes negative, leading to the instability of soliton. By taking 
min(Im(cj 2 2{±±})) = 0, one obtains 

_ ay/ 5a 2 - AA^J^ 2 -a 2 - 4 7 2 5a 2 - AA^j 2 -a 2 - 4 7 2 ^ 
£cr ~ 2(/l+ VT 2- ^) 

which yields an approximate boundary for the onset of instability, e.g., in the 
(e, a)-plane for fixed A and 7. 



2{A+J^ 2 



3.2.2 Intersite type II 

From our analysis of Eqs. (12) and (13), we obtain the discrete eigenvalues 
for an intersite bright soliton type II, which are given, with errors of order e 2 , 
by 



= i \a ± y 4/1 \J~i 2 - a 2 - A^ 2 + 5a 2 ± 



\jAA\fif~- a 2 - Aj 2 + 5a 2 / 
(47) 



^ 2 ( TT | =i\a± lA/lyV -a 2 - 4 7 2 + 5a 2 ± 2{ A ^2 ] 
\ Y AAyJ-f 2 - a 2 - Aj 2 + 5a 2 J 

(48) 

assuming the term (AA^/j 2 — a 2 — A-f 2 + 5a 2 ) is 0(1). Notice that wi{ T=F } 
and W2{zpip} are equal when a = \J A^ 2 — A 2 /2. 
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When a < 7, both min(Im(u;i and min(Im(o;2{ =F=F })) are negative at 
£ = and always decrease as e is increased; the decrement of min(Im(cj2{ TT })) 
is greater than min(Im(a;i{ =F=F })) for a > \/4:j 2 — A 2 /2. When a = 7, 
min(Im(cji^ TT i.)) and mm(lm(uj2{ TT })) arc zero at e = 0. At nonzero e, 
the former remains zero, but the latter becomes negative and decreases as e 
increases. These facts allow us to conclude that an intersite bright soliton type 
II is always unstable, except at a = 7 and e = 0. One can check that when 
a = 7, the eigenvalues of intersite type II are the same as in intersite type I. 



3.2.3 Intersite type III and IV 

As intersite type III and IV are symmetric, their eigenvalues are exactly the 
same. Our calculation shows the following. 

(i) For the case a < a t h, the eigenvalues of the intersite type III and IV, up 
to C(e 2 ), are 



w h{±t} = w ii{t±} = *a ± V 4/1^2 - a 2 + A-f 2 - ha 2 

(2 7 ^7 2 - a 2 - Aj + ay 'A 2 - 7 2 + a 2 )e 

± / 1 

7^ 4/1^2 - a 2 + 47 2 - 5a 2 

w i2{± T } = w i2{ T ±} = « ± \ f±A V7 2 -"a 2 " _ 4 7 2 + 5 a 2 

± (2 7 ^7 2 - a 2 + yl 7 - qyg2 _ 7 2 + a 2) £ ^ ^ ^ 

7^ 4/1 V7 2 - a 2 - 47 2 + 5a 2 
(ii) For the case a t h < ce < 7, the eigenvalues, up to order e 2 , arc 



w 2i{± T } = ^2i{ T ±} = i (a ± ^/ 5a 2 - 4/1^ 7 2 - a 2 - 47 

t (27/^2—^2- _ Al + Qv/yl 2 _ 7 2 + a 2) £ N ( _ | ( 

5a 2 - AAy/y* - - a 2 - 4"f~ 



^22{± T } = ^22{=|=±} =i[a±\J 4A\/"f 2 - a 2 - 4-y 2 + 5a 2 

± (27^72 - a 2 + Ay - a^J A 2 - 7 2 + a 2 )e \ 
AA^/y 2 - a 2 - 47 2 + 5a 2 / 

We should assume again that the terms ±(AAy/j 2 — a 2 + 4j 2 — 5a 2 ) and 
(4/1 y^ 2 — a 2 — A'y 2 + 5a 2 ) in the above expansions are of 0(1). 
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In the first case, the eigenvalues (50) are apparently pure imaginary, with 
an imaginary part whose minimum value is negative for all e. In the second 
case, it is clear that for a < 7 the minimum value of the imaginary part 
of the eigenvalues (51) is positive (less than a) initially at e — and then 
increases as e increases. However, for this case (a < 7), the minimum value 
of the imaginary part of the eigenvalues (52), which are exactly the same as 
the eigenvalues (50), is negative at e = and then decreases as e is turned 
on. In contrast, for a — 7 the minimum value of the imaginary part of the 
eigenvalues (51) and (52) remains zero for all e. The above fact shows that 
both intersite soliton type III and IV are always unstable, except at a = 7. 
In fact, as shown in the numerical calculation later, the intersite type III and 
IV are no longer defined along this line, due to a pitchfork bifurcation with 
intersite type I. 

4 Comparisons with numerical results, and bifurcations 

In order to find the numerical solutions for each soliton discussed in the 
previous section, we solve the stationary equation (2) [cf. Eq. (4)] using 
a Newton-Raphson (NR) method. The evaluation is performed in domain 
n e [—N,N], i.e., for a lattice of 27V + 1 sites, with periodic boundary con- 
ditions u±(jv+i) = u tN- As an initial guess, we use the corresponding exact 
soliton solutions in the uncoupled limit e = from which we then numerically 
continue for nonzero e. As an illustrative example, the numerical solutions 
for each type of onsite and intersite bright soliton with parameter values 
(e, yl,7, a) = (0.1,1,0.5,0.1) are depicted in Fig. 1. The corresponding ana- 
lytical approximations are also plotted therein showing good agreement with 
the numerical results. 

To examine the stability of each soliton, we solve the eigenvalue problem 
(6) numerically and then compare the results with the analytical calculations. 
Moreover, we show later that the relevant solitons experience saddle-node 
and/or pitchfork bifurcations. To depict the diagram of these bifurcations, we 
use a pseudo-arclength method which allows us to continue the solution past 
turning points (by varying one parameter). In addition, our analysis of the 
eigenvalues for some particular solutions leads to the fact of the presence of 
Hopf bifurcations. We will determine the nature of Hopf bifurcation points 
and perform continuation of the bifurcating limit cycles in the next section 
by employing the numerical continuation package Matcont. 

In all illustrative examples below, we use N — 50 which is large enough to 
capture the behavior of the soliton in an infinite domain but not too costly in 
numerical computations. In addition, for the sake of simplicity, we set A = 1 
and 7 = 0.5. 
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(a) Onsite type I (b) Onsite type II 




n n 

(c) Intersite type I (d) Intersite type II 




(e) Intersite type III (f) Intersite type IV 



Fig. 1. Profiles of onsite and intersite bright solitons of different types, as indicated 
in the caption of each panel, for parameter values (e, A, 7, a) — (0.1, 1, 0.5, 0.1). Solid 
lines show the numerical results while dashed lines indicate the analytical approxi- 
mations given by Eqs. (22) and (23) for the onsite type I and II, respectively, and 
by Eqs. (34), (35), (36), and (37) for the intersite type I, II, III, and IV, respec- 
tively. The circle and cross markers correspond to the real and imaginary part of 
the solutions, respectively. 
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4.1 Onsite bright solitons 
4.1.1 Onsite type I 

We start by testing the validity of our analytical approximation for the critical 
eigenvalues given by Eqs. (29) and (30). We present in Fig. 2 comparisons 
between the analytical and numerical results for the critical eigenvalues as 
functions of e. We plot comparisons for three values a = 0.1,0.47,0.497 to 
represent the cases a < a cp , a cp < a < a t h, and a t h < a < 7, respectively 
(see again the relevant discussion in the previous section). From the figure, 
we conclude that our prediction for small e is relatively close to the numerics. 

2 1 , , , , 1 

1.8 - ^ 

1.6 - 
1.4 - 
_ 1 .2 

3_ 

1 - 

0.8 - — 

0.6 - 

0.2 1 ' ' ' ' 

0.02 0.04 0.06 0.08 0.1 

e 

Fig. 2. Comparisons between the critical eigenvalues of an onsite bright soliton 
type I obtained numerically (solid lines) and analytically (dashed lines). The upper 
and middle curves correspond, respectively, to a = 0.1 and a = 0.485, which are 
approximated by Eq. (29), whereas the lower corresponds to a = 0.497, which is 
approximated by Eq. (30). 

For the three values of a given above, we now present in Fig. 3 the eigen- 
value structure of the soliton and the corresponding diagram for the imaginary 
part of the critical eigenvalues as functions of e. Let us now describe the results 
in more detail. 

First, we notice that at e = the critical eigenvalues for a = 0.1 lie be- 
yond the outer band of the continuous spectrum, while for a = 0.485 they 
are trapped between the two inner bands of the continuous spectrum. As s is 
turned on, the corresponding critical eigenvalues for a = 0.1 and a = 0.485 
collide with, respectively, the outer and the inner bands, leading to the bi- 
furcation of the corresponding eigenvalues. The minimum imaginary part of 
these bifurcating eigenvalues, however, does not immediately become negative. 
Hence, for relatively small e we conclude that the soliton is always stable; this 
in accordance with our analytical prediction of the previous section. The crit- 
ical values of e at which min(Im(w)) = indicating the onset of the instability 
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are depicted by the star markers in panels (c) and (f) in Fig. 3. Interestingly, 
for a — 0.485 there is a re-stabilization of the soliton as shown by the larger 
£ star marker in panel (f). 

Next, for a = 0.497 the discrete eigenvalues initially (at s = 0) lie on the 
imaginary axis; they come in pairs and are symmetric about the line Im(cj) = 
a = 0.497, furthermore the minimum one is above the real axis. When e 
increases, both eigenvalues approach one another and finally collide at the 
point (0, a = 0.497) creating a new pair of discrete eigenvalues along the line 
Im(w) = a = 0.497. Each pair of the eigenvalues then again bifurcates after 
hitting the inner edge of the continuous spectrum. However, the minimum 
imaginary part of these bifurcating eigenvalues is always greater than zero 
even for larger e [sec panel (i)]. From this fact, we therefore conclude that the 
soliton in this case is always stable. This conclusion agrees with our analytical 
investigation. 

The minimum value of Im(w) (in color representation) of the onsite bright 
soliton type I for a relatively large range of e and a gives the (in)stability 
region in the (e, a)-plane as presented in Fig. 4. The stable region is indeed 
determined whenever min (Im(w)) > for each e and a. The lower and upper 
dotted horizontal lines in this figure, i.e., respectively, a — a cp rts 0.4583 and 
a = ctth ~ 0.49659, represent the boundaries of the regions which distinguish 
the description of the eigenvalue structure of the soliton. The solid line in this 
figure indicates the (in)stability boundary, i.e., when min(Im(w)) = 0. Three 
representative points (star markers) lying on this line reconfirm the corre- 
sponding points in panels (c) and (f) in Fig. 3. As shown in the figure, there 
is an interval of a in which the soliton is stable for all e. This is interesting as 
the onsite soliton, which was shown [11, 12] to be destabilized by a parametric 
driving, now can be re-stabilized by a damping constant. 

Let us revisit Fig. 3 for a = 0.1 and a = 0.485. We notice that at zero- 
crossing points e c (shown by the star markers in panels (c) and (f)), the 
following conditions hold: 

(i) There is a pair (equal and opposite) of non-zero real eigenvalues, and 

(ii) The e-derivative of the imaginary part of the pair of eigenvalues mentioned 
in (i) is non-zero at e c . 

The second condition is also called the transversality condition. We assume 
that the so-called first Lyapunov coefficient of the zero-crossing points is 
nonzero, i.e. the genericity condition. According to the Hopf bifurcation the- 
orem (see, e.g., Ref. [29], keeping in mind that our eigenvalue is denoted by 
iuj), the above conditions imply that at e = e c Eq. (1) has time-periodic (limit 
cycle) solutions bifurcating from a (steady-state) onsite bright soliton type I. 
We then call such a critical point e c a Hopf point. By applying the centre 
manifold theorem, for example, we can generally determine the nature of a 
Hopf point e c through its first Lyapunov coefficient h(e c ) (see, e.g., Ref. [29]); 
the Hopf bifurcation is subcritical iff h(e c ) > and supercritical iff l\{e c ) < 0. 
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Fig. 3. The first and second columns of panels show the (Re(aj), Im(aj))-plane of the 
eigenvalues of onsite bright solitons type I for several values of a and e, as indicated 
in the caption of each panel (each row of panels depicts three different values of a). 
For a — 0.1 and a = 0.485, the corresponding left and middle panels illustrate the 
eigenvalues of stable and unstable solitons. The third column shows the path of the 
imaginary part of the critical eigenvalues uj clit as functions of e for the corresponding 
a. The locations of e at which Im(w C rit) = are indicated by the star markers. 



Because the occurrence of Hopf bifurcation in the onsite type I also in- 
dicates the onset of (in)stability, the collection of Hopf bifurcation points in 
the (e, a)-plane therefore lies precisely on the (in)stability boundary line (see 
again Fig. 4). However, at the stationary point e w 1.46 the condition (ii) 
for the occurrence of a (non-degenerate) Hopf bifurcation does not hold. At 
this special point, we have a saddle-node bifurcation of Hopf points, i.e. a 
double-Hopf (Hopf-Hopf) bifurcation. Due to the violation of the transver- 
sality condition, there may be no periodic solution or even multiple periodic 
solutions at the denegerate point. We will examine this point later in Sec. 5, 
where it will be shown through numerical continuations of limit cycles near 
the degenerate point that the former possibility occurs. 
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Fig. 4. (Color online) The (in)stability region of onsite bright solitons type I in the 
(e, a)-plane. The corresponding color represents the minimum value of Im(cj) (for all 
eigenvalues ui) for each e and a. Thus, the region in which min (Im(aj)) > indicates 
the region of stable soliton, otherwise unstable. The boundary of stable-unstable 
regions, i.e., when min(Im(oj)) = 0, is given by the solid line (three representative 
points (star markers) on this line correspond to those points in panels (c) and (f) 
in Fig. 3). The boundary curve also indicates the occurrence of Hopf bifurcations 
with one degenerate point, i.e. a double-Hopf bifurcation, at e ~ 1.46 as indicated 
by the white-filled circle. The lower and upper horizontal dotted lines correspond to 
Eqs. (31) and (26), respectively (see text). 



4.1.2 Onsite type II 

For this type of solution, a comparison between the critical eigenvalues ob- 
tained by analytical calculation, which is given by Eq. (32), and by numerics, 
is presented in Fig. 5. We conclude that our analytical prediction for small e 
is quite accurate. 

The eigenvalue structure of onsite solitons type II for a = 0.1 and the two 
values e = 0.1, 1 and the corresponding curve of imaginary part of the critical 
eigenvalues are given in Fig. 6. This figure shows that the soliton is always 
unstable even for a large e. This fact is consistent with the analytical predic- 
tion. We notice in the figure that there is a new pair of discrete eigenvalues 
bifurcating from the inner edge of continuous spectrum at relatively large e 
[see panel (b)]. 

By evaluating the minimum value of Im(cj) for a relatively large e and a, we 
obtain that the soliton is always unstable for a < 7 = 0.5 and, contrastingly, 
stable for a — 7. In the latter case, the eigenvalues of the onsite type II are 
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3 0.915 




0.002 0.004 0.006 0.008 0.01 



Fig. 5. Comparison between the critical eigenvalues of onsite bright solitons type II 
for a — 0.1 produced by numerics (solid line) and by analytical approximation (32) 
(dashed line). 
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Fig. 6. The top panels show the eigenvalue structure of onsite bright solitons type 
II for a = 0.1 and two values of e as indicated in the caption. The bottom panel 
depicts the imaginary part of the critical eigenvalues as a function of e. 
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exactly the same as in the onsite type I; the minimum value of the imaginary 
part remains zero for all e. 

4.1.3 Saddle-node bifurcation of onsite bright solitons 

We observed from numerics and analytics that when approaching a = 7, the 
onsite bright soliton type I and type II possess the same profile as well as 
the same stability, consistent with the saddle-node bifurcation experienced 
by the two solitons. A diagram of this bifurcation can be produced, e.g., by 
plotting the norm of the numerical solution of these two solitons as a function 
of a for fixed s = 0.1. To do so, we apply a pseudo-arc-length method to 
perform the numerical continuation, starting from the onsite type I at a = 0. 
The obtained diagram is presented in Fig. 7 and the corresponding analytical 
approximation is also depicted therein. As shown in the figure, the onsite 
type I, which is stable, turns into the onsite type II, which is unstable. Both 
numerics and analytics give the same turning point [or so-called limit point 
(LP)] at a = 7 = 0.5. We also conclude that the analytical approximation for 
the norm is quite close to the numerics, with the accuracy for the onsite type 
I better than type II. Indeed, their accuracy could be improved if one uses 
smaller e. 




Fig. 7. A saddle-node bifurcation of onsite bright solitons for e = 0.1. The onsite 
type I (stable) merges with the onsite type II (unstable) at a limit point (LP) a = 
7 = 0.5. The solid and dashed lines represent the norm of the solutions obtained by 
numerical calculation and analytical approximation, respectively. The insets depict 
the profile of the corresponding solutions at the two values a — 0.1, 0.5. 
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4.2 Intersite bright solitons 
4.2.1 Intersite type I 

Let us first compare our analytical prediction for the critical eigenvalues, 
given by Eqs. (42)-(43) and (44)-(45), with the corresponding numerical re- 
sults. We present the comparisons in Fig. 8 by considering three values of 
a = 0.1,0.465,0.497 as representative points for the three cases discussed in 
the previous section. From the figure we see that the double eigenvalues which 
coincide originally at e — then split into two distinct eigenvalues as e in- 
creases. We conclude that our approximation for small e is generally quite 
accurate. 




0.15 
0.1 
0.05 



0.005 0.01 0.015 0.02 

e 

(c) a = 0.497 

Fig. 8. Comparisons of the two distinct critical eigenvalues of intersite bright soli- 
tons type I obtained numerically (solid lines) and analytically (dashed lines) for 
three values of a as indicated in the caption of each panel. The upper and lower 
curves in panels (a) and (b) are plotted from, respectively, Eqs. (42) and (43), while 
in panel (c) from Eqs. (45) and (44). 



Next, we move on to the description of the eigenvalue structure of the 
intersite bright solitons type I and the corresponding imaginary part of the 
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two critical eigenvalues as functions of e; these are depicted in Fig. 9 for the 
three values of a used before. The first and second columns in the figure rep- 
resent conditions of stability and instability, respectively. For a = 0.1, the two 
critical eigenvalues successively collide with the outer band of the continuous 
spectrum and the corresponding bifurcating eigenvalues coming from the first 
collision contribute to the instability. For a — 0.465, the first collision is be- 
tween one of the critical eigenvalues with the inner edge of the continuous 
spectrum. The second collision is between the other critical eigenvalue with 
its pair. In contrast to the previous case, the instability in this case is caused 
by the bifurcating eigenvalues coming from the second collision. Moreover, 
for a = 0.497, contribution to the instability is given by one of the critical 
eigenvalues moving down along the imaginary axis. All the numerical results 
described above are in accordance with our analytical observations in Sec. 3. 

Let us now focus our attention on the right panels of Fig. 9 by particularly 
discussing the properties of the critical points of e at which the curve of the 
minimum imaginary part of the critical eigenvalues crosses the real axis (these 
are shown by the star markers). The first and third points (from left to right) in 
panel (c) as well as the points in panels (f ) and (i) indicate the onset of stable- 
to-unstable transition. Contrastingly, the second point in panel (c) illustrates 
the beginning of the re-stabilization of solitons. In fact, the first three points 
in panel (c) mentioned above admit all conditions for the occurrence of a Hopf 
bifurcation (see again the relevant explanation about these conditions in our 
discussion of onsite type I); therefore, they also correspond to Hopf points. In 
addition, the fourth point of zero crossing in panel (c), which comes from one 
of the purely imaginary eigenvalues, indicates the branch point of a pitchfork 
bifurcation experienced by the solutions of intersite type I, III, and IV. We 
will discuss this type of bifurcation in more detail in the next section. 

The (in)stability region of intersite bright solitons type I in the (e, em- 
plane is given by Fig. 10. In the figure, we also depict the two distinguish- 
able (solid and dashed) lines representing the two distinct critical eigenvalues 
whose imaginary parts become zero. The star points on the lines correspond 
to those points in the right panels of Fig. 9. The boundary line which sepa- 
rates the stable and unstable regions in the figure is shown by the bold (solid 
and dashed) lines. The lower and upper dotted horizontal lines in the figure 
represent, respectively, a = a cp w 0.4583 and a — ath ~ 0.49659 which di- 
vide the region into three different descriptions of the eigenvalue structure. 
Interestingly, for a t h < a, we can make an approximation for the numerically 
obtained stability boundary (see the inset). This approximation is given by 
Eq. (46) which is quite close to the numerics for small e. 

We notice in Fig. 10 that the solid line (not the rightmost) and dashed 
line also represent Hopf bifurcations, with one special point (the white-filled 
circle) which does not meet the second condition for the occurrence of a (non- 
degenerate) Hopf bifurcation mentioned above. We will analyse the special 
point in the next section. We see from the figure that the bold parts of the 
Hopf lines coincide with the (in) stability boundary, while the nonbold ones 
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Fig. 9. The first and second columns of panels show the structure of eigenvalues 
of intersite bright solitons type I in the complex plane, for three values of a, each 
of which uses two different values of e, to depict the condition of stability (left 
panel) and instability (middle panel). The third column shows the imaginary part 
of the two distinct critical eigenvalues as functions of £ for the corresponding a. The 
locations of zero-crossings in these panels are indicated by the star markers. 



exist in the unstable region. In addition, we also observe that the rightmost 
solid line in Fig. 10 indicates the collection of branch points of pitchfork 
bifurcation experienced by the intersite type I, III, and IV; the bold part of 
this line also indicates the (in) stability boundary. 

4.2.2 Intersite type II 

For intersite bright solitons type II, we present in Fig. 11 a comparison of 
two critical eigenvalues between the numerics and the analytical calculation 
given by Eqs. (47) and (48). We see from the figure that our approximation 
for relatively small s is quite close to the numerics. The snapshot of the eigen- 
value structure of this type of solution for two points (a, e) and the path of 
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Fig. 10. (Color online) As Fig. 4 but for intersite bright solitons type I. The bound- 
ary between stable and unstable regions is given by the bold (solid and dashed) 
lines. The dashed-dotted line in the inset is our analytical approximation given by 
Eq. (46). The Hopf bifurcation lines are depicted by the solid (not the rightmost) and 
dashed lines. The white-filled circle indicates a degenerate Hopf point. The branch 
points of pitchfork bifurcation are shown by the rightmost solid lines. 

the imaginary part of corresponding two discrete eigenvalues are depicted in 
Fig. 12. We conclude that the intersite soliton type II is unstable even for 
large e. 

Moreover, the evaluation of the minimum value of Im(w) of the intersite 
bright solitons type II in the (e, a)-plane gives the (in)stability window (not 
shown here). It is shown that the soliton, except at the point a = 7 = 0.5 and 
£ = 0, is always unstable. This result agrees with our analytical prediction. 

4.2.3 Intersite type III and IV 

Now we examine the intersite bright soliton type III which, due to symmetry, 
has exactly the same eigenvalues as type IV. Shown in Fig. 13 is the analytical 
approximation for two critical eigenvalues given by Eqs. (49)-(50) or (51)-(52), 
which are compared with the corresponding numerical results. We conclude 
that the approximation is quite accurate for small e and that the range of 
accuracy is wider for smaller value of a. 

The structure of the eigenvalues of this type of solution and the curves 
of the imaginary part of the corresponding two critical eigenvalues are given 
by Fig. 14 for the three values of a used in Fig. 13. The figure reveals the 
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Fig. 11. Comparison of two critical eigenvalues of intersite bright solitons type II for 
a = 0.3 between numerics (solid lines) and analytics (dashed lines). The analytical 
approximation is given by Eq. (47) (lower curve) and Eq. (48) (upper curve). 
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Fig. 12. (a)-(b) The eigenvalue structure of intersite bright solitons type II for two 
values (a, e) as indicated in the caption, (c) The imaginary part of two discrete 
eigenvalues in varied e. 
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condition of instability of solitons up to the limit points of e at which the 
minimum imaginary part of the eigenvalues becomes zero; these conditions are 
indicated by the corresponding vertical lines in the third column. In fact, these 
limit points indicate the branch points of pitchfork bifurcation experienced by 
the intersite solitons type I, III, and IV (we will discuss this bifurcation in 
more detail in the next section). 

The first and second columns of Fig. 14 respectively present the condition 
just before and after a collision of one of the discrete eigenvalues which does 
not contribute to the instability of solitons. Interestingly, as shown in panel 
(c) , such an eigenvalue also crosses the real axis at some critical e as indicated 
by the empty circle. The latter condition, in fact, indicates a Hopf bifurcation, 
which occurs when the soliton is already in unstable mode. This is different 
from the previous discussions where the Hopf bifurcations also indicate the 
change of stability of solitons. 



0.005 0.01 0.015 0.02 0.025 0.03 



0.002 0.004 0.006 0.008 0.01 



(a) a = 0.1 



(b) a = 0.465 




(c) a = 0.497 



Fig. 13. Comparisons between two critical eigenvalues of intersite bright solitons III 
and IV obtained numerically (solid lines) and analytically (dashed lines) for values 
of a as shown in the caption. In panels (a) and (b), the upper and lower dashed 
curves correspond, respectively, to Eqs. (49) and (50), whereas in panel (c) they 
correspond to Eqs. (51) and (52). 
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Presented in Fig. 15 is the (in)stability window for intersite bright solitons 
type III and IV which is defined as the area to the left of the solid line; this 
line represents the set of the branch points of pitchfork bifurcation in the 
(e, a)-plane. From the figure, we conclude that the intersite type III and IV 
are always unstable. The area to the right of the solid line belongs to the 
unstable region of intersite type I. One can check that this line is exactly 
the same as the rightmost solid line in Fig. 10. In addition, the dashed line 
appearing in Fig. 15 depicts the occurrence of Hopf bifurcations. However, 
there is one special point indicated by the white-filled circle, at which the 
e-derivative of the imaginary part of the corresponding critical eigenvalue is 
zero; this degenerate point will be discussed further in Sec. 5. The empty 
circle lying on the Hopf line reconfirms the corresponding point in panel (c) 
of Fig. 14. 

4.2.4 Saddle-node and pitchfork bifurcation of intersite bright 
solitons 

From both numerical and analytical results discussed above, we observed that 
the intersite type I and type II have the same profile and stability when ap- 
proaching a = 7. This fact indicates the appearance of a saddle- node bifur- 
cation undergone by the two solitons. Moreover, there also exists a pitchfork 
bifurcation experienced by the intersite type I, III, and IV. 

One can check that the norm of the intersite type III and IV is exactly 
the same for all parameter values so that this quantity can no longer be used 
for depicting a clear bifurcation diagram. Therefore, we now simply plot the 
value of |uo| 2 for each solution, e.g., as a function of a and fixed e = 0.1; this 
is shown in Fig. 16 where the numerics (solid lines) is obtained by a pseudo- 
arc-length method. As seen in the figure, the intersite type I, III, and IV meet 
at a (pitchfork) branch point (BP) a ss 0.49. At this point, the stability of 
the intersite type I is switched. Furthermore, the intersite type I and II also 
experience a saddle-node bifurcation where they merge at a limit point (LP) 
a = 7 = 0.5. Just before this point, the intersite type I possesses one unstable 
eigenvalue, while the type II has two unstable eigenvalues. The two critical 
eigenvalues for the intersite type I and II then coincide at LP. We confirm 
that our analytical approximation for the value of |iio| 2 is relatively close to 
the corresponding numerical counterpart. 

Next, let us plot the value of |uo| 2 for each soliton by fixing a = 0.1 and 
varying e (presented in Fig. 17). The pitchfork bifurcation experienced by the 
intersite type I (solid line), type III (upper dashed line), and type IV (lower 
dashed line) is clearly shown in the figure. The three solitons meet together 
at a branch point BP. We also depict in the figure the points at which Hopf 
bifurcations emerge (labelled by indexed H). For the shake of completeness, 
we also plot the relevant curve for the intersite type II (dotted line). 
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Fig. 14. (First and second columns) The structure of eigenvalues of intersite bright 
solitons type III and IV for parameter values (a, e) as indicated in the caption. 
(Third column) The imaginary part of two critical eigenvalues obtained by varying 
e. The vertical lines indicate the limit points of e up to which the soliton exists, i.e., 
when the minimum imaginary part of the eigenvalues becomes zero. 



5 Nature of Hopf bifurcations and continuation of limit 
cycles 

If there is only one pair of non-zero real eigenvalues and the other eigenval- 
ues have strictly positive imaginary parts, a Hopf bifurcation also indicates 
the change of stability of the steady state solution. In this case, the periodic 
solutions bifurcating from the Hopf point coexist with cither the stable or un- 
stable mode of the steady state solution. If the periodic solutions coexist with 
the unstable steady state solution, they are stable and the Hopf bifurcation is 
called supercritical. On the other hand, if the periodic solutions coexist with 
the stable steady state solution, they are unstable and the Hopf bifurcation 
is called subcritical. 
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Fig. 15. The (in)stability region of intersite bright solitons type III and IV in 
(e,a)-space. The solid line indicates the branch-point line of pitchfork bifurcation. 
The dashed line represents the occurrence of Hopf bifurcations (with one degenerate 
point at the white-filled circle), which arises from one of the critical eigenvalues 
which does not contribute to the instability of solitons. The empty circle lying on 
the dashed line corresponds to that point depicted in panel (c) of Fig. 14. 

To numerically calculate the first Lyapunov coefficient for a Hopf point 
and perform a continuation of the bifurcating limit cycle, we use the numerical 
continuation package Matcont. Due to the limitations of Matcont, we evaluate 
the soliton using 21 sites. In fact, this setting does not affect significantly the 
soliton behavior compared to that used in the previous section. 

In this section, we examine the nature of Hopf points and the stability of 
cycle continuations in onsite type I, intersite type I, and intersite type III-IV. 

5.1 Onsite type I 

For this type of solution, in particular at a = 0.1, we have one Hopf point, 
which occurs at e c « 0.3077 (see again panel (c) in Fig. 3). From Matcont, 
we obtain l\(e c ~ 0.3077) > which indicates that the Hopf point e c is sub- 
critical and hence the limit cycle bifurcating from this point is unstable. A 
continuation of the corresponding limit cycle is given in Fig. 18(a). As the 
Hopf point in this case also indicates the change of stability of the station- 
ary soliton, one can confirm that the bifurcating periodic solitons are stable 
because they coexist with the stable onsite type I; this agrees with the com- 
puted first Lyapunov coefficient above. Interestingly, the continuation of the 
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Fig. 16. Saddle-node and pitchfork bifurcations of intersite bright solitons by vary- 
ing a and fixing e = 0.1. The curves depict the value of |ito | 2 of each solutions 
obtained numerically (solid lines) and analytically (dashed lines). The profiles of 
the corresponding solutions at some values of a are shown in the relevant insets. 
The intersite type I, III, and IV merge at a branch point (BP) a ~ 0.49 and the 
intersite type I and II meet at a limit point (LP) a — 7 = 0.5. 

limit cycle also experiences saddle-node and torus bifurcations, as indicated 
by the points labelled LPC (limit point cycle) and NS (Neimark-Sacker), re- 
spectively. The profile of a representative periodic soliton over one period is 
shown in Fig. 18(b), from which we clearly see the typical oscillation in the 
soliton amplitude. 

From the previous discussion we have mentioned that there is one degen- 
erate point for Hopf bifurcations in onsite type I, which is indicated by the 
white-filled circle in Fig. 4. In Fig. 19, we depict numerical continuations of pe- 
riodic orbits of two Hopf bifurcations near the degenerate point. We obtained 
that the limit cycle branches bifurcating from the Hopf points are connected 
and form a closed loop. This informs us that as a approaches the critical value 
for a degenerate Hopf point, the "radius" of the loop tends to zero. Hence, 
one may conclude that at the doublc-Hopf point, there is no bifurcation of 
periodic orbits. 

5.2 Intersite type I 

In particular for a = 0.1, there are three Hopf points detected for the inter- 
site type I (see again Fig. 17). For point Hi (e w 0.2782), Matcont gives a 
negative value for the first Lyapunov coefficient, which means that the bifur- 
cating periodic soliton is stable or Hi is supercritical. The corresponding cycle 
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Fig. 17. A pitchfork bifurcation of intersite bright solitons for fixed a = 0.1 and 
varied e. The curves represent the numerical value of |tto| 2 for the corresponding 
solutions as a function of e. The intersite type I (dashed line), type III (upper 
dashed line), and type IV (lower dashed line) merge at a branch point (BP). The 
occurrence of Hopf bifurcation (Hi) is detected in intersite type I, III, and IV. The 
dotted line corresponds to the intersite type II. 




0.25 0.26 0.27 0.28 0.29 0.3 0.31 0.32 

e 



(a) (b) 



Fig. 18. (a) The continuation of the limit cycle from a Hopf point H for an onsitc 
soliton type I with a — 0.1. The first Lyapunov coefficient for H calculated by 
Matcont is positive, i.e. H is subcritical. The bold solid line represents the norm 
of the stationary soliton while the solid and dashed lines indicate, respectively, the 
maximum and minimum of the norm of the bifurcating periodic solitons. (b) The 
profile of an unstable periodic soliton (as H is subscritical) over one period (T ~ 
3.4319) corresponding to the black-filled circle in panel (a). 
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Fig. 19. As Fig. 18(a) but for a = 0.492642. Two Hopf points (stars) in the neigh- 
bourhood of the degenerate point (the while-filled circle in Fig. 4) are shown to be 
connected by a branch of limit cycles. 



continuation is presented in Fig. 20(a). As shown in the figure, the limit cy- 
cle bifurcating from Hi coexist with the unstable mode of the (steady-state) 
intersite type I which confirms the supercritical Hi. This is valid because 
the Hopf bifurcation in this case also indicates the change of stability of the 
soliton. We also see from the figure that the cycle continuation contains NS, 
LPC, and BPC (branch point cycle) points which indicate the occurrence 
of, respectively, torus, saddle-node, and pitchfork bifurcations for limit cycle. 
The branches of the cycle continuation from the BPC point are shown in the 
figure. A representative periodic soliton (in one period) which occurs at one 
representative point along the cycle continuation is depicted in Fig. 20(b), 
which shows the oscillation between the two excited sites. 

Next, for H 2 (e w 0.3871) and H 3 (e w 0.4934), the first Lyapunov coeffi- 
cients given by Matcont are negative and positive valued, respectively. Thus, 
H 2 is supercritical while H 3 is subcritical, which implies that the limit cycle 
bifurcating from H 2 and H 3 are stable and unstable, respectively. The contin- 
uations of the corresponding limit cycles are shown in Fig. 21(a). From the 
figure, we see that the limit cycles bifurcating from H 2 and H 3 respectively 
coexist with the unstable and stable stationary intersite soliton type I. This 
fact is consistent with the nature of H 2 and H 3 as given by Matcont. In ad- 
dition, as shown in the figure, a period-doubling (PD) bifurcation also occurs 
in the cycle continuation coming from H 3 . This bifurcation seems to coincide 
with the turning point of cycle (LPC) which appears in the cycle continua- 
tion starting from H 2 . The profile of one-period periodic solitons at the two 
representative points near H 2 and H 3 are presented in Figs. 21(b) and 21(c), 
respectively. We cannot see clearly the typical oscillation of the periodic soli- 
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Fig. 20. (a) The cycle continuation from Hopf point Hi for intersite bright soliton 
type I with a — 0.1. In this case, Hi is supercritical. The bold solid line indicates the 
value of |wo| 2 for the stationary soliton, which is the same as that shown in Fig. 17. 
The solid and dashed lines represent, respectively, the maximum and minimum value 
of |«o| 2 for the bifurcating periodic solitons, which also experience a pitchfork cycle 
bifurcation. The branches of the cycle are depicted by the dash-dotted (maximum 
|ito| 2 ) and dotted (minimum |wo| 2 ) lines, (b) The profile of a stable periodic soliton 
(as Hi is supercritical) over one period (T w 5.5265) corresponding to the star point 
in panel (a) . (c,d) Enlargements of, respectively, the upper and the lower rectangles 
in panel (a). 



ton in Fig. 21(b) as it occurs very near to H 2 . By contrast, the oscillation in 
the soliton amplitude is clearly visible in Fig. 21(c). 

Similarly to the onsite type I, we also noticed the presence of a double- 
Hopf bifurcation in the intersite type I, i.e., the white-filled circle in Fig. 10. 
To investigate the point, we evaluate several Hopf points nearby the bifur- 
cation point and perform numerical continuations for limit cycles, which are 
presented in Fig. 22. Unlike the case in the onsite type I, here the (non- 
degenerate) Hopf points are not connected to each other by a closed loop of a 
branch of limit cycles. As we observe this scenario at any Hopf point that is 
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(b) (c) 

Fig. 21. (a) As Fig. 20(a) but for H2 and H3, where the inset gives the zoom-in for 
the corresponding region showing that H2 and H3 are supercritical and subcritical, 
respectively. The bold solid line is the same as that shown in Fig. 17, i.e., representing 
the value of |wo| 2 for the stationary intersite soliton. The solid (dashed) and dash- 
dotted (dotted) lines shows the maximum (minimum) value of |uo| for the periodic 
soliton which bifurcates from, respectively, H2 and H3. (b,c) The profile of periodic 
solitons over one period T « 6.708 and T « 3.4985 which corresponds, respectively, 
to the star and the black-filled circle in panel (a). From the nature of H2 and H3, 
periodic solitons in (b) and (c) are stable and unstable, respectively. 
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arbitrarily close (up to a numerical accuracy) to the degenerate (codimension 
2) bifurcation, it indicates that at the double-Hopf point, there is a bifurcation 
of at least two branches of periodic solutions. 




Fig. 22. As Fig. 21(a) but for a = 0.108 (triangles) and a = 0.11082 (stars) in the 
proximity of the while-filled circle in Fig. 15. 



5.3 Intersite type III and IV 

As intersite bright soliton type III and IV possess the same eigenvalue struc- 
tures, the nature of the corresponding Hopf bifurcation and the stability of 
the continuation of each limit cycle will be the same as well. Therefore it is 
sufficient to devote our discussion to intersite type III only. 

As shown in Fig. 17, there is one Hopf point, namely H4, for the intersite 
type III at a = 0.1. In this type of solution, the Hopf bifurcation occurs while 
other eigenvalues already give rise to instability; this is different from the 
type of Hopf bifurcation discussed previously. Therefore we cannot perform 
the analysis as before in determining the stability of the bifurcating periodic 
soliton. In fact, according to calculation given by Matcont, the first Lyapunov 
coefficient for H4 is positive (subcritical), which means that the bifurcating 
periodic soliton is unstable. 

Fig. 23(a) shows the continuation of the corresponding limit cycle from 
H4. A representative one-period periodic soliton at e near H4 (indicated by 
the black- filled circle) is shown in Fig. 23(b), from which we can see clearly 
the oscillation in the amplitude of soliton. 

Next, we study the double- Hopf bifurcation for the intersite type III- IV 
shown by the white-filled circle in Fig. 15. Presented in Fig. 24 is the contin- 
uation of the limit cycles from two Hopf points about the degenerate point, 
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Fig. 23. (a) The cycle continuation from Hopf point H4 for intersite bright soliton 
type III with a = 0.1 showing that H4 is subcritical. The dashed line shows the 
value of |mo| for the stationary soliton (the same as that shown in Fig. 17) while 
the solid and dotted lines represent, respectively, the maximum and minimum |mo| 2 
of the bifurcating periodic solitons. (b) The profile of an unstable periodic soliton 
over one period T « 3.7388 corresponding to the black-filled circle in panel (a). 



from which we see that they are connected to each other. Therefore, as for 
the case of the onsite type I, we argue that there is no bifurcation of periodic 
solutions at the degenerate point. 




Fig. 24. As Fig. 23(a) but for a = 0.1411. The degenerate (Hopf- Hopf) point 
indicated as the while-filled circle in Fig. 15 is at e ~ 0.49. 
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6 Conclusion 

In this paper, we have considered a parametrically driven damped discrete 
nonlinear Schrodingcr (PDDNLS) equation. The existence and stability of 
fundamental discrete bright solitons have been examined analytically through 
a perturbation theory for small e and then corroborated by numerical calcu- 
lations. We showed that there are two types of onsite discrete soliton, namely 
onsite type I and II. For onsite type I, we found an interval in a for which the 
soliton is stable for any coupling constant, i.e., a damping can rc-stabilize a 
driven onsite soliton. Contrastingly, the onsite type II was found to be always 
unstable for all e. These two solitons experience a saddle- node bifurcation 
with the limit point a — 7 for any e. 

We also showed that there are four types of intersite discrete soliton, called 
intersite type I, II, III, and IV. In fact, intersite type III and IV are essentially 
considered as one solution due to its symmetry. We obtained that intersite type 
I in the region of instability in the non-dissipative case can be stabilized by 
damping while intersite type II and III-IV are always unstable. A saddlc-nodc 
bifurcation, as for the onsite soliton, was found to be undergone by intersite 
type I and II. Moreover, we also obtained that intersite type I, III, and IV 
experience a pitchfork bifurcation. The branch points of such a bifurcation in 
the (e, a)-plane have been calculated numerically. 

More interestingly, we observed that Hopf bifurcation also occurs in onsite 
type I, intersite type I, and intersite type III-IV, which confirms the existence 
of the corresponding periodic solitons (limit cycles) in the PDDNLS equation. 
The continuation of the limit cycles as well as the stability of the periodic 
solitons have been demonstrated numerically using the numerical continuation 
software Matcont. In particular, subcritical Hopf bifurcations for onsite type 
I and intersite type III-IV were observed. Moreover, we obtained three Hopf 
bifurcations for intersite type I. It was shown that two of these points generate 
stable periodic solitons, i.e. the bifurcations are supercritical. 

Note that similar studies for the continuum limit of Eq. (1) have been put 
forward in, e.g., Refs. [24, 25, 26]. Hopf bifurcations and the corresponding 
periodic solitons were reported and discussed therein. The connection between 
the results presented in this work, which correspond to weakly coupled lattices, 
and those of [24, 25, 26] are proposed for future study. 
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